Noiseless scattering states in a chaotic cavity 
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Shot noise in a chaotic cavity (Lyapunov exponent A, level spacing 5, linear dimension L), coupled 
by two TV-mode point contacts to electron reservoirs, is studied as a measure of the crossover from 
stochastic quantum transport to deterministic classical transport. The transition proceeds through 
the formation of fully transmitted or reflected scattering states, which we construct explicitly. The 
fully transmitted states contribute to the mean current 7, but not to the shot-noise power S. We flnd 
that these noiseless transmission channels do not exist for N < ^JkpL, where we expect the random- 
matrix result 5/26/ = 1/4. ForiV > y/kpL we predict a suppression of the noise oc (kpL/N^)^^^^'^'^ . 
This nonlinear contact dependence of the noise could help to distinguish ballistic chaotic scattering 
from random impurity scattering in quantum transport. 

PACS numbers: 05.45.Mt, 03.65.Sq, 72.70.+m, 73.63.Kv 



Shot noise can distinguish deterministic scattering, 
characteristic for particles, from stochastic scattering, 
characteristic for waves. Particle dynamics is determin- 
istic: A given initial position and momentum fixes the 
entire trajectory. In particular, it fixes whether the par- 
ticle will be transmitted or reflected, so the scattering 
is noiseless. Wave dynamics is stochastic: The quan- 
tum uncertainty in position and momentum introduces a 
probabiUstic element into the dynamics, so it is noisy. 

The suppression of shot noise in a conductor with de- 
terministic scattering was predicted many years ago from 
this qualitative argument 1]. A better understanding, 
and a quantitative description, of how shot noise mea- 
sures the transition from particle to wave dynamics in a 
chaotic quantum dot was put forward by Agam, Aleiner, 
and Larkin 0, and developed further in Ref. 0. The 
key concept is the Ehrenfest time te, which is the char- 
acteristic time scale of quantum chaos The noise 
power S cx cxp^—te/td) was predicted to vanish expo- 
nentially with the ratio of te and the mean dwell time 
td = nh/NS in the quantum dot (with 6 the level spac- 
ing and N the number of modes in each of the two point 
contacts through which the current is passed). A recent 
measurement of the A^-dependence of S is consistent with 
this prediction for te < to, although an alternative ex- 
planation in terms of short-range impurity scattering de- 
scribes the data equally well 5]. 

The theory of Ref. introduces the stochastic element 
by means of long-range impurity scattering and adjusts 
the scattering rate so as to mimic the effect of a finite 
Ehrenfest time. Here we take the alternative approach 
of explicitly constructing noiseless channels in a chaotic 
quantum dot. These are scattering states which are ei- 
ther fully transmitted or fully reflected in the semiclas- 
sical limit. They are not described by random matrix 
theory By determining what fraction of the avail- 
able channels is noiseless, we can deduce a precise upper 
bound for the shot-noise power. A random matrix con- 
jecture for the remaining noisy channels gives an explicit 
form of S{N). We find that the onset of the classical sup- 




FIG. 1: Selected equipotentials of the electron billiard. The 
outer equipotential is at Ep, the other equipotentials are at 
increments of 0.19 Ep. Dashed lines a and h show the sections 
described in the text. Also shown is a flux tube of transmitted 
trajectories, all originating from a single closed contour in 
a transmission band, representing the spatial extension of a 
fully transmitted scattering state. The flux tube is wide at 
the two openings and squeezed inside the billiard. 

pression of the noise is described not only by the Ehren- 
fest time, but by the difference of r^; and the ergodic time 
To, which we introduce and calculate in the paper. The 
resulting nonlinear dependence of In 5* on A'^ may help to 
distinguish between the competing explanations of the 
experimental data |^. 

We illustrate the construction of noiseless scattering 
states for the two-dimensional billiard with smooth con- 
fining potential U{x,y) shown in Fig. 1. The outer 
equipotential defines the area in the x — y plane which is 
classically accessible at the Fermi energy Ee ~ p\/2m 
(with Pf = T^kE the Fermi momentum). The motion in 
the closed billiard is chaotic with Lyapunov exponent A. 
We assume the billiard to be connected at a; = and 
X = Lhy two similar point contacts to leads of width W 
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FIG. 2: Section of phase space at px = \/p'p — Py and x — 0, 
corresponding to line a in Fig. 1. Each dot in this surface 
of section is the starting point of a classical trajectory that 
is transmitted through the lead at a; = L (black/red), or 
reflected back through a; = (gray/green). The points lie 
in narrow bands. Only the trajectories with dwell time t < 
12mL/pF are shown. 

extended along the ix-direction. 

The beam of electrons injected through a point con- 
tact into the billiard has cross-section W and transverse 
momenta in the range {—pw,Pw)- The number of chan- 
nels N ~ pwW/h in the lead is much smaller than the 
number of channels M ~ ppL/h supported by a typical 
cross-section of the billiard. While W/L ^ 1 in general, 
the ratio pw /pf depends on details of the potential near 
the point contact, li pw/PF ^ 1 one speaks of a col- 
limated beam. This is typical for a smooth potential, 
while a hard- wall potential typically has pw pf (no 
collimation) . We define r,„iii = miTi{W/L,p\Y /pf) and 
^max = m.a.x{W/L,pw/pF)- 

The classical phase space is four-dimensional. By re- 
stricting the energy to Ep and taking a; = we obtain 
the two-dimensional section of phase space shown in Fig. 
121 The accessible values of y and Py lie in a disc-shaped 
region of area A = Nh in this surface of section. Up to 
factors of order unity, the disk has width rmin and length 
''max (if coordinate and momentum are measured in units 
of L and pp, respectively). In Fig. 2 one has rmin — ''max- 
Each point in the disc defines a classical trajectory that 
enters the billiard (for positive Px) and then leaves the bil- 
liard either through the same lead (reflection) or through 
the other lead (transmission). The points lie in narrow 
bands, which we will refer to as "transmission bands" 
and "reflection bands" . 

It is evident from Fig. |21 that the area Aj enclosed by 
a typical transmission (or reflection) band j is much less 
than A. For an estimate we consider the time t{y,py) 
that elapses before transmission. Let tj be the dwell time 
averaged over the starting points y^Py in a single band. 
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FIG. 3: Dwell-time distribution for the billiard of Fig. 1. 
Electrons at the Fermi energy are injected through the left 
lead. Time is in units of mL/pp- Inset: the same data on 
a semilogarithmic scale with larger bin size of the histogram. 
Three characteristic time scales are seen: tw, tq, td- 

The fluctuations of t around the average are of the order 
of the time tw — mW/pw to cross the point contact, 
which is typically <C tj. As we will see below, the area 
of the band decreases with tj as 

~ ^oexp(-Atj) if tj :$> 1/X,tw- (1) 

The prefactor = ^''min/''max depends on the degree of 
collimation. In Ref. the symmetric case rmin = '"max 
was assumed, when = 

We now proceed to the construction of fully transmit- 
ted scattering states. To this end we consider a closed 
contour C within a transmission band j. The starting 
points on the contour define a family of trajectories that 
form a flux tube inside the billiard (see Fig. 1). The 
semiclassical wave function 

fp(x,y) = ^ \/ pa{x,y) exp[iSa{x,y)/h] (2) 

is determined as usual from the action So- and density pa- 
that solve the Hamilton- Jacobi and continuity equations 

\VSf = 2m{Ep -U), V • (pV5) = 0. (3) 

The action is multivalued and the index a labels the dif- 
ferent sheets. Typically, there are two sheets, one origi- 
nating from the upper half of the contour C and one from 
the lower half. 

The requirement that ijj is single-valued as one winds 
around the contour imposes a quantization condition on 
the enclosed area, 

j Pydy = {n + l/2)h. (4) 

The increment 1/2 accounts for the phase shift acquired 
at the two turning points on the contour. The quantum 
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number n = 0, 1, 2, . . . is the channel index. The largest 
value of n occurs for a contour enclosing an area Aj. 
The number of transmission channels Nj within band j 
is therefore given by Aj/h, with an accuracy of order 
unity. In view of Eq. |^ we have 



Nj ~ {Ao/h)exp{-Xtj), for tj 



N, = 0, 



for 



tj > TE- 



(5a) 
(5b) 



The time 



TE = \n{Ao/h) = ln(A^r„,i„/r,„ax) (6) 

above which there are no fully transmitted channels is 
the Ehrenfest time of this problem. 

By decomposing one of these Nj scattering states into 
a given basis of transverse modes in the lead one con- 
structs an eigenvector of the transmission matrix product 
tt^ . The corresponding eigenvalue 7),„ equals unity with 
exponential accuracy in the semiclassical limit n 3> 1. 
Because of the degeneracy of this eigenvalue any linear 
combination of eigenvectors is again an eigenvector. This 
manifests itself in our construction as an arbitrariness in 
the choice of C . 

We observe in Fig. 1 that the spatial density profile 
p{x,y) of a fully transmitted scattering state is highly 
non-uniform. The flux tube is broad (width of order W) 
at the two openings, but is squeezed down to very small 
width inside the billiard. A similar effect was noted 
in the excited states of an Andreev billiard (a cavity con- 
nected to a superconductor). Following the same argu- 
ment we estimate the minimal width of the flux tube as 



Ly/N-j/kpL. 

The total number 



=Y,N, = N r p{t)dt 

Jo 



(7) 



of fully transmitted and reflected channels is determined 
by the dwell-time distribution P{t) 0- Fig. 3 shows this 
distribution in our billiard. One sees three different time 
scales. The narrow peaks represent individual transmis- 
sion (reflection) bands. They consist of an abrupt jump 
followed by an exponential decay with time constant tw- 
These exponential tails correspond to the borders of the 
bands, where the trajectory bounces many times between 
the sides of the point contact. If we smooth P{t) over 
such short time intervals, an exponential decay with time 
constant td = nh/NS is obtained (inset). The decay 
starts at the so called "ergodic time" tq. There are no tra- 
jectories leaving the cavity for t < To. So the smoothed 
dwell-time distribution has the form. 



P{t) = CXp[(To - t)/TD]eit - To), 



(8) 



with 9{t) the unit step function. 

In order to flnd tq we consider Fig. 4, where the sec- 
tion of phase space along a cut through the middle of the 
billiard is shown (line b in Fig. 1). It is convenient to mea- 
sure the momentum and coordinate along b in units otpp 




FIG. 4: Section of phase space in the middle of the billiard, 
along line b in Fig. 1. The subscript || indicates the component 
of coordinate and momentum along this line. Elongated black 
areas Oj show the positions of the 5-th crossing of the injected 
beam with this surface of section. The area Oinitiai is the 
position of the first crossing. Points inside Ofinai leave the 
billiard without further crossing of line b. For times less than 
the ergodic time tq there is no intersection between Oj and 

Oflnal • 

and L. The injected beam crosses the section for the first 
time over an area Oinitiai of size Tmax x ''min = hN/ppL. 
(Fig. 4 has r^in — ^'max, but the estimates hold for any 
Tmin < ''max < !■) Further crossings consist of increas- 
ingly more elongated areas. The fifth crossing is shown 
in Fig. 4. The flux tube intersects line 6 in a few disjunct 
areas Oj, of width r-mine^'*'* and total length rmaxe^*. 
(Due to conservation of the integral ^ p • dr enclosing 
the flux tube, the total area J2j Oj decreases only when 
particles leave the billiard.) The typical separation of ad- 
jacent areas is (rmaxC'^*)"^. To leave the billiard (through 
the right contact) without further crossing of 6 a particle 
should pass through an area Ofinai — ''max x ''min- This 
is highly improbable Q until the separation of the areas 
Oj becomes of order rmax, leading to the ergodic time 



To 



A^Mnr: 



The ergodic time varies from to ^ A ^ 

To = \~'^\n{kFL/N) for r^in ^ ''max- 



(9) 

for Tniax ^ 1 tO 

The overlap of 



the areas Oj and Ofinai is the mapping of the transmis- 
sion band onto the surface of section b. It has an area 



PpLrl 



-At 



A{rn 



1 /''max)e leading to Eq. p. 
Substituting Eq. JHJ) into Eq. |(7J) we arrive at the num- 
ber TVq of fully transmitted and reflected channels. 



No = N9{te - To) 



I _ g(TO-TE)/TB 



(10) 

TE-To = X-^lniN^/kFL). (11) 
There are no fully transmitted or reflected channels if 
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te < To, hence if TV < ^/kpL. Notice that the depen- 
dence of Te and tq separately on the degree of coUima- 
tion drops out of the difference te — tq. The number of 
noiseless channels is therefore insensitive to details of the 
confining potential. An Ehrenfest time oc \n{N'^ /kpL) 
has a ppe ared before in connection with the Andrccv bil- 
liard 10], but the role of coUimation (and the associated 
finite ergodic time) was not considered there. 

Eqs. m and ((SJ imply that the majority of noiseless 
channels group in bands having Nj ^ 1, which justifies 
the semiclassical approximation. The total number of 
these noiseless bands is {N — Nq) / Xtd , which is much less 
than both N — Nq and Nq. Because of this inequality the 
relatively short trajectories contributing to the noiseless 
channels are well separated in phase space from other, 
longer trajectories (cf. Fig. 2). 

The shot-noi se p ower S is related to the transmission 
eigenvalues by |iy 

N 

S = 2eIg''Y.'^k{l-%), (12) 

fc=i 

with / the time- averaged current and 5 = X^A; 
dimcnsionless conductance. The Nq fully transmitted or 
reflected channels have 7^, = 1 or 0, hence they do not 
contribute to the noise. The remaining N — Nq channels 
contribute at most 1/4 per channel to Sg/2el. Using 
that g = N/2 for large iV, we arrive at an upper bound 
for the noise power S < e/(l — Nq/N). 

For a more quantitative description of the noise power 
we need to know the distribution P{T) of the transmis- 
sion eigenvalues for the N — Nq noisy channels, which 
can not be described semiclassically. We expect the 
distribution to have the same bimodal form P{T) = 
7r-ir-i/2(l - r)-i/2 as in the case iVo = This ex- 
pectation is motivated by the earlier observation that the 
iVo noiseless channels are well separated in phase space 
from the N — Nq noisy ones. Using this form of P{T) we 



find that the contribution to Sg/2el per noisy channel 
equals JqT{1 - T)P{T)dT = 1/8, half the maximum 
value. The Fano factor F — S/2eI is thus estimated as 

F = \, iov N < ^/k^, (13a) 
F = {{kpL/N^)'^^'''^'^, for N>y/k^. (13b) 

This result should be compared with that of 
Ref. 0: F' = likpL)''^^/''^^ . The ratio F' / F = 
eyiY>[{2N5/TTh\)h\{N/kpL)\ is always close to unity (be- 
cause NS/nhX ~ N/kpL < 1). But F - j and 
F' — ^ are entirely different for N ^ VkpL, which 
is the relevant regime in the experiment 0. There 
the A^-dependence of the shot noise was fitted as F = 
j{l — tQ/rjj) — |(1 — constant x A^), where tg is some 
A^-independent time. Eq. p3|l predicts a more com- 
plex A^-dependence, a plateau followed by a decrease as 
InF (X —Nln{N'^/kFL), which could be observable if the 
experiment extends over a larger range of A^. 

We mention two other experimentally observable fea- 
tures of the theory presented here. The reduction of the 
Fano factor described by Eq. is the cumulative effect 
of many noiseless bands. The appearance of new bands 
with increasing A^ introduces a fine structure in F(N), 
consisting of a series of cusps with a square-root singu- 
larity near the cusp. The second feature is the highly 
nonuniform spatial extension of open channels, evident 
in Fig. 1, which could be observed with the STM tech- 
nique of Ref. . From a more general perspective the 
noiseless channels constructed in this paper show that 
the random matrix approach may be used in ballistic sys- 
tems only for sufhciently small openings: A^ ^ y/kpL is 
required. For larger N the scattering becomes determin- 
istic, rather than stochastic, and random matrix theory 
starts to break down. 
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